Epitaxial interfaces between crystallographically mismatched materials 
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We report an unexpected mechanism by which an epitaxial interface can form between materials having 
strongly mismatched lattice constants. A simple model is proposed in which one material tilts out of the interface 
plane to create a coincidence-site lattice that balances two competing geometrical criteria — low residual strain 
and short coincidence-lattice period. We apply this model, along with complementary first-principles total- 
energy calculations, to the interface formed by molecular-beam epitaxy of cubic Fe on hexagonal GaN and find 
excellent agreement between theory and experiment. 

PACS numbers: 81.15.Aa, 81.15.Hi, 68.55.-a, 68.37.Lp 



A fundamental goal of materials science is to eluci- 
date and exploit the physical principles that govern epi- 
taxial growth 111 01 ■ Some of these principles are well- 
established. For example, if the lattice constants of the film 
and substrate are close but not identical then a coherently 
strained film may grow up to a critical thickness, beyond 
which misfit dislocations relieve the strain [3|]. Alterna- 
tively, a film and substrate having lattice spacings close to 
an integer ratio m/n may form an epitaxial interface de- 
scribed by a coincidence lattice Etl • 

In this Letter we report a new and unexpected mecha- 
nism by which epitaxial films can grow on substrates hav- 
ing, in principle, an arbitrary lattice mismatch. We illus- 
trate this mechanism experimentally by growing single- 
crystal Fe on M-plane GaN. The Fe grows in an unusual 
crystallographic orientation with a very high Miller index, 
Fe(205). We develop a simple theoretical model which, 
when complemented with total-energy calculations, cor- 
rectly predicts this exact orientation as well as the single- 
domain nature of the film. Finally, we use our model to 
propose a new strategy for growing nonpolar epitaxial GaN 
films on high-index Si substrates. 

The epitaxial growth of both GaN and Fe was performed 
in a custom-built molecular-beam epitaxy system equipped 
with solid-source effusion cells for Ga and Fe. Active 
nitrogen was provided by a radio-frequency N2 plasma 
source. Nucleation and growth were monitored in situ 
by reflection high-energy electron diffraction. A 130-nm 
thick layer of M-plane GaN was first grown on a 4H- 
SiC(l 100) substrate under Ga-stable conditions and a tem- 
perature of 720°C. After growth of the GaN layer, excess 
Ga was desorbed prior to cooling down to 350°C for the 
deposition of Fe |6]. The Fe film grew at this temperature 
at a rate of 0.13 nm/min to a final thickness of 27 nm. The 
resulting Fe/GaN/SiC heterostructure was investigated by 
cross-sectional transmission electron microscopy (TEM) 
and convergent-beam electron diffraction (CBED) using a 
JEOL JEM-3010 operating at 300 keV. Electron backscat- 




FIG. 1: (a) Cross-sectional transmission electron micrograph 
taken along the [TT20] zone axis of the Fe/GaN/4H-SiC structure 
under investigation. The top 2 nm of the Fe film are oxidized, 
(b) Cross-sectional high-resolution transmission electron micro- 
graph of the Fe/GaN(l 100) interface along the [TT20] zone axis. 
The Fe(101) and GaN(0001) lattice planes form an angle, a, of 
approximately 23°. Stacking faults in the GaN, visible in both 
panels, do not affect the orientation of the Fe film. 



tering diffraction (EBSD) was carried out in a Zeiss Ultra- 
55 scanning electron microscope equipped with an EDAX- 
TSL EBSD system. 

Figure [Tf a) shows a cross-sectional transmission elec- 
tron micrograph of the Fe/GaN/4H-SiC structure. Despite 
the structural and chemical differences of the constituent 
materials, the respective layers are well defined and ex- 
hibit abrupt interfaces. The GaN layer is seen to contain 
stacking faults due to the stacking mismatch between 2H- 
GaN and 4H-SiC. Nevertheless, the high-resolution detail 
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FIG. 2: Schematic view of a commensurate interface between 
a film with lattice constant df and a substrate with lattice con- 
stant d s . The interface plane is (h0£) with respect to the film, 
corresponding to a tilting from (001) toward (101) by the angle 
6 = tan - (h/£). The condition for commensurability is that a 
coincidence-site lattice, defined by a pair of integers (m,n), ex- 
ists such that md s = n(h 2 + £ 2 ) l l 2 df. The example depicted here 
is (205) with 6 = 21.8° (equivalent to the angle a = 23.2° as 
defined in Fig.fj} and (m,n) = (3, 1). 
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FIG. 3: Misfit strain in commensurate interfaces between bcc Fe 
and GaN(l 100). The points represent all possible ¥e(h0£) planes 
having h < I < 75. The misfit e xx is the strain component along 
the [0001] direction of the GaN. The size of each plot symbol 
is inversely proportional to the period m of the coincidence-site 
lattice (CSL) that minimizes the misfit. Fe(MM) planes having 
both small misfit strain and small CSL period are labeled. 



of the Fe/GaN interfacial region shown in Fig.QJb) reveals 
an epitaxial relationship between the Fe film and the under- 
lying GaN(lTOO) layer. The (T01) lattice planes of the Fe 
film are clearly resolved and are found to be well ordered, 
unaffected by the stacking disorder in the GaN layer. Of 
special interest is the angle, approximately 23°, formed by 
the Fe (101) planes and the vertical interface normal. This 
angle indicates that the Fe interface plane has a high Miller 
index — an unexpected finding in light of the comparatively 
high surface energies of high-index metal surfaces. We 
show now that precisely this orientation is predicted by a 
simple, physically transparent model (complemented with 
first-principles total-energy calculations) of epitaxial inter- 
faces between dissimilar materials. 

Consider the formation of an interface between a film 
(/) and a substrate (s) having different lattice constants df 
and d s . If the lattice mismatch is sufficiently small then 
the strained film may grow coherently until it reaches its 
critical thickness fl. For much larger mismatch this sce- 
nario becomes unlikely. Epitaxial growth is nevertheless 
possible by tilting the orientation of the film, as the TEM 
image in Fig. |TJb) makes clear. Figure [2] illustrates how 
an arbitrary lattice mismatch can be accommodated by al- 
lowing the film to have an orientation between (001) and 
(101) given by the Miller indices (h0£). Our goal below is 
to predict the most stable film orientation given the lattice 
constants df and d s . 

The unit cell of a film with orientation (h0£) has length 
Lf = (h 2 + £ 2 ) x l 2 df. In order for the film and substrate 
to be commensurate there must exist a coincidence-site 



lattice (CSL), defined by a pair of integers (m,n), such 
that md s = nLf. This condition is unrealistically restric- 
tive, however. In real systems the film will tolerate a small 
compressive or tensile strain e xv which relaxes the CSL 
condition to md s = nLf{ \ + e VT ). 

We propose two geometrical criteria for identifying can- 
didate orientations for interfaces with low energy. (1) The 
misfit strain £ xx should be as small as possible, and (2) the 
period m of the CSL should be as small as possible. The 
latter criterion is motivated by analogy to low-energy grain 
boundaries between two identical materials, which often 
have a CSL with small unit cell volume £ J3l. For inter- 
faces between different materials it is not generally possi- 
ble to minimize the strain and CSL period simultaneously. 
Nor is it obvious how to construct a single objective func- 
tion of both which could then be optimized. Instead, we 
apply both criteria with the aim of selecting a small sub- 
set of candidate orientations for subsequent study with a 
more quantitative method such as density-functional the- 
ory (DFT). 

To apply these criteria to the growth of Fe on the M- 
plane of GaN we equate df with the bcc Fe lattice constant, 
2.866 A 0], and d s with the GaN c lattice parameter, 5.186 
A fl]. Figure [3] shows the resulting Fe misfit strain needed 
to satisfy the CSL condition for a large number of hypo- 
thetical orientations of the Fe film. In this plot each ori- 
entation (/i0£) is represented by its angle 9 — tan~ l (h/l) 
relative to the (001) plane. The period of each CSL is en- 
coded by the size of the plot symbol, which is inversely 
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FIG. 4: (a) Pseudo-Kikuchi pattern of the Fe film on GaN(l 100) 
depicted in Fig.Q] The major low-index zone axes are indicated, 
(b) In-plane orientation map of the Fe film. Variations in color 
encode local deviations from the nominal orientation. The film 
has the same orientation over the entire area within a tolerance of 
1°. (c) Model of the Fe(205)/GaN(lT00) epitaxial interface. For 
clarity an ideal geometry with arbitrary registry is shown here. 
Slightly less than one unit cell along [0001] is depicted. A con- 
ventional unit cell of Fe is highlighted in yellow. 



proportional to m. Only points with small strains, less than 
3%, are displayed here. 

The rich structure visible in Fig. [3] makes it clear that 
the misfit strain can be made arbitrarily small for many 
different film orientations. Hence the strain alone cannot 
provide a definitive criterion favoring a particular growth 
plane. Moreover, the vast majority of these low-strain ori- 
entations require a very large CSL period and hence do not 
constitute physically meaningful commensurability. Only 
a very few orientations offer both a small strain and small 
CSL period, namely (207), (205), and (203). The (205) 
orientation corresponds to the angle = 21.8°. This is 
equivalent to a = 45 — = 23.2° as defined in Fig. [TJb) 
and thus is in excellent agreement with the measured an- 
gle, 23°, obtained from TEM. 

It is important to realize that our purely geometrical cri- 
teria do not distinguish between the orientation (205) de- 



picted in Fig. |2] and its symmetry-equivalent counterpart, 
(205) (0 = -21.8°), created by rotating the film by 180° 
about the substrate normal. Indeed, the pairs [(h0£), (hOl)] 
have the same strain and CSL and hence are equivalent 
within this model. In the Fe/GaN system, however, the 
polarity of the GaN wurtzite structure breaks this equiva- 
lence. The question that then arises is whether the influ- 
ence of the polarity is sufficiently strong to select a sin- 
gle orientation, and if so, which one? The TEM image 
in Fig. Q] shows a single orientation but is limited to a 
nanometer-scale region of the film. To characterize a much 
larger area we used EBSD |9J]. Figure Ufa) shows the re- 
sulting pseudo-Kikuchi pattern of the Fe film in Fig.Q] The 
pattern exhibits sharp and well defined Kikuchi bands, re- 
flecting high crystal quality of the Fe film and allowing for 
a fast and reliable indexing of the patterns recorded while 
scanning the electron beam over a large area (9 Aim 2 ). The 
resulting EBSD map shown in Fig. @fb) visualizes the in- 
plane orientation of the Fe film with a spatial and angular 
resolution of 20 nm and 1°, respectively. The map reveals 
the complete absence of any domain structure. Indeed, the 
film is single crystalline, and has the same orientation as 
found by TEM throughout the mapped area. 

To understand why a single orientational domain is 
found requires going beyond a model based solely on in- 
terface geometry. Now the interface structure — its precise 
atomic arrangement and chemical bonding — must be ad- 
dressed. To do this we used DFT to calculate the relative 
formation energies of finite Fe films grown on the M-plane 
of GaN. We considered the six different orientations pre- 
dicted by the geometrical model to be favorable: (207), 
(205), (203), and their rotated counterparts (207), (205), 
(203). The different films were the same thickness, 6 A, 
equivalent to about four monolayers. Figure Uc) depicts 
the Fe(205)/GaN(l 100) interface as an example. 

We have previously shown that at 3-4 monolayers the 
contribution of the Fe/GaN interface formation energy to 
the full formation energy of the film is already converged 
We also find that the Fe free surface energy varies 
by less than 1 meV/A 2 among the three orientations we 
consider here ill ill . Therefore the relative formation energy 
of the finite film closely mirrors, with good accuracy, the 
formation energy of the isolated interface. 

For each orientation the Fe film was slightly strained 
along the x [GaN(0001)] direction according to its CSL 
as discussed above. There is also a lattice mismatch in 
the y [GaN(1120))] direction because the Fe lattice con- 
stant and GaN a lattice parameter differ by nearly 12%. 
This mismatch was accommodated by a single CSL, com- 
mon to all orientations, containing eight unit cells of Fe 
and seven of GaN. The GaN substrate was represented 
by a slab of four atomic layers with fixed in-plane equi- 
librium lattice parameters and a passivating bottom layer. 
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TABLE I: Formation energies calculated within density- 
functional theory of Fe/GaN(1100) films for the Fe planes in 
Fig. |3]that have both a small misfit strain and a coincidence-site 
lattice (CSL) with small period m. The Fe film thickness is the 
same for all cases, approximately 6 A. Formation energies are 
relative to the most favorable plane, Fe(205), in units of meV/A 2 . 

Plane Misfit strain CSL period Formation energy 
Fe(205) +0.008 3 

Fe(205) +0.008 3 4 

Fe(203) +0.004 2 9 

Fe(203) +0.004 2 16 



polar GaN films. The basic idea is to turn Fig. [2] upside- 
down and consider the growth of M-plane GaN on a suit- 
able high-index substrate. For a given material this re- 
quires identifying candidate orientations corresponding to 
small strain and small CSL period. One very promising 
material is Si, which is already in widespread use as a 
flat substrate for GaN/Si epitaxy despite the resulting high 
dislocation densities lfl6ll . Many high-index Si substrates 
are readily available, and some have already been used for 
GaN growths 117111811 . Calculations are in progress to iden- 
tify promising high-index Si orientations for growing non- 
polar GaN with low strain fl^l . 
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Total energies and forces were calculated within the PBE 
generalized-gradient approximation Ill2ll to DFT using 
projector-augmented- wave potentials as implemented in 



VASP A 1 3L 1 1411 . All Fe and GaN atomic positions were 
fully relaxed except the bottom GaN layer. For each Fe 
film orientation the formation energy depends strongly on 
the choice of Fe-GaN interface registry. We systematically 
varied the registry over a grid in both x and y to locate the 
global energy minimum for each orientation. The plane- 
wave cutoff for all calculations was 400 eV. 

The resulting formation energies are listed in Table U 
The most favorable orientation is Fe(205). Of the candi- 
dates tested, this orientation has neither the smallest possi- 
ble strain nor the shortest possible CSL period, indicating 
that the optimal interface structure is an important third 
criterion that must supplement the two geometrical crite- 
ria. Note also that the small variation in the formation en- 
ergy of the Fe surface, which is included in the formation 
energy of the film, is too small to affect the overall energy 
ordering of the orientations. 

To distinguish Fe(205) from Fe(205) experimentally re- 
quires determining the absolute polarity of the GaN sub- 
strate. We did this by recording CBED patterns with a 
beam spot size of approximately 15 nm under two-beam 
conditions. Simulations of CBED patterns were performed 
using JEMS II 1 511 to index the crystallographic directions 
observed in the experimental patterns and thus to deter- 
mine the polarity. The resulting absolute orientation rela- 
tionship is Fe(205)[0T0]||GaN(lT00)[TT20], in agreement 
with the prediction of our theoretical model complemented 
with the results of DFT calculations. 

Our model suggests a new strategy for growing non- 
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